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Abstract 



The lattice technique of studying the strong interaction of matter is used to obtain predictions 
of the hadronic spectrum. These simulations were performed by the UKQCD collaboration 
using full (unquenched) QCD. Details of the results, a comparison with quenched data, and 
novel methods of extracting spectral properties are described. 



1 Introduction 



The lattice technique of studying field theories such as QCD (the strong force of particle physics) 
has had a long and successful history [|J. In this method, the usual computational modelling trick of 
discretising space-time, and replacing derivatives by finite differences is employed. This allows the 
study of field theories at any value of the coupling, and contrasts with perturbation theory which is 
restricted to systems where the interaction strength is weak. For this reason, the lattice technique 
is particularly applicable to calculations involving the strong interaction of particle physics, and is 
now regarded as the method of choice for its study. It has the advantage over of its competitors in 
that it has no model assumptions, and is also systematically improvable. 

Historically, due to computational resource issues, lattice simulations have used the "quenched" 
approximation (which ignores virtual quark loops in the vacuum). Now that computer power is in 
the Teraflop region, lattice QCD studies are able to simulate without this approximation. 

The UKQCD collaboration^ has played a leading role in the study of QCD via the lattice tech- 
nique. Most of the present emphasis of this collaboration is the study of hadrons using "full" (i.e. 
"unquenched" ) QCD. A summary of recent results will be given (sec [||, [J for details) and an 
exciting new approach to obtaining spectral information from the lattice based on the "Maximum 
Entropy Method" will be discussed. 

The main computer resource of the UKQCD collaboration is currently a Cray T3E in Edinburgh, 
however this will soon be strengthened by the acquisition of an APEmille computer in Swansea. 
Furthermore, funding has been approved for a supercomputer in the multi- Teraflop range in the 
medium term. 

2 Brief Overview of Lattice Gauge Theory 

In lattice QCD, the quark degrees of freedom are defined on the sites of the 4-dimensional lattice, 
and the gluonic degrees of freedom (corresponding to the carriers of the strong force) are defined 
on the links. This enables a lattice version of the theory to be defined which maintains as many of 
the symmetries of the continuum theory as possible. (For a reviews of lattice QCD see ||.) 

Typical quantities of interest can then be expressed as n-point functions, G n (t), of hadronic oper- 
ators. Specifically these are weighted averages over configurations of the degrees of freedom with 
a Boltzmann-like weight, e~@ s where S is the lattice version of the continuum QCD action. It is 
important to note that this is exactly analogous to statistical mechanics. This means that many 
of the techniques developed in that area are applied to lattice QCD - for example the use of the 
Monte Carlo integration method. 

Quantities of physical interest, such as hadronic masses, can be extracted from calculations of the 
2-point function G%{t) as follows. It can be shown that 

G 2 (t) =Y J Z i e~ m -> Z e- Mot as t -> oo. (1) 

i 

The sum is over all states \i > with the same quantum numbers as the operators which define the 
2-point function, G2(t). Hence the ground state mass of a hadron, Mq, can be extracted simply by 
fitting G2{t) to an exponential at sufficiently large times t. 

This is not the end of the story though because Mq (and any other quantity extracted from a 
lattice calculation) is a function of the input parameters of the simulation. These are the coupling 
strength go, the quark mass m and the volume L 3 . Ideally what is required is the continuum limit 
(a — > 0) of Mq which corresponds to go — > (since QCD is asymptotically free). However, as has 
been found by several groups (e.g. §| H|), the quark mass m also has impact on the lattice spacing 
a, i.e. a = f(go,m). This will have an important implication as we will see in a moment. 
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Another issue regarding lattice QCD is the error which is introduced in the discretisation of the 
continuum theory. In the conventional approach, the lattice action reproduces the continuum action 
only up to terms of 0(a). Over the last several years, lattice groups, including UKQCD, have been 
using so-called improved versions of the lattice action which have lattice errors of 0(a 2 ). While 
this is an obvious improvement over conventional lattice actions, it does not preclude lattice results 
being contaminated with lattice artefacts, all-be-it at the 0(a 2 ) level. 

In order to establish a clean signal for many studies it is imperative to fix a = const for two reasons: 
(i) as seen above, the lattice calculation leaves the answer correct to 0(a 2 ), so in order to keep 
this error fixed, a must obviously be fixed; (ii) in order to maintain fixed finite-size effects, a must 
be fixed (simply because the lattice size L = Na where N is the number of lattice sites along a 
side). Thus, considering that a = f(go,m) (see above) the requirement that a = const defines a 
trajectory in the (go,m) plane. 

3 Results 

In the previous section we argued that simulating along a fixed a trajectory in (go,m) parameter 
space lead to the cleanist study. UKQCD have followed this philosophy and performed simulations 
at four different parameter choices for (go,m) such that a = const ~ 0.1 fm (see j|] for details). Note 
that one of these four parameter choices corresponds to the quenched approximation (specifically 
m = oo). It is important to note that when a physical quantity changes value between these four 
parameter values, the cause can be attributed to variations in the quark mass m (i.e. "unquenching" 
effects), and not lattice systematics such as finite volume effects (see previous section). 

In this section an overview of three quantities of particular interest is given. 
3.1 Static Quark Potential 

One of the fundamental quantities that can be studied on the lattice is V(r), the potential between 
two infinitely heavy quarks (the "static quark potential"). Phenomenologically, 

V{r) = V + - + ar. (2) 
r 

In Figure 1, V(r) is plotted for all four parameter choices. The phenomenological model Eq.(||) is 
also shown. As can be seen there is agreement amongst the data and Eq.(^) to a high precision. 
However, close inspection of the data near the origin shows a systematic tendency towards a higher 
value of a as the quark mass m is decreased towards it physical value. We take this as a signal for 
unquenching effects [f§ . 
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Figure 1: Static quark potential versus distance r on left, and Hadronic masses My versus Mp s 
on right. All quantities have been scaled by the distance scale rg. 



3.2 Hadronic Spectrum 

In this subsection the hadronic spectrum is discussed. We concentrate here on the pseudoscalar 
and vector mesonic states whose masses are denoted My and Mps respectively, since these can be 
determined most unambiguously. Theoretically we expect the following relationship to hold: 

My = A + B Ml s , 

where A and B are constants. In Figure 1, My is plotted against Mp S for all four parameter 
choices. Note that within each parameter choice, the "valence" quark mass is still a free parameter 
hence the multiple data points for each parameter choice in the figure Q. Also plotted in Figure 1 
is the experimental values for the ir, p, K, K* , rj s and </> mesons @. As can be seen there is a trend 
towards the experimental values as the quenched approximation is removed. 



3.3 Topological Susceptibility 



The topological susceptibility, x> plotted in Figure 2, shows a marked unquenching effect. This 
is interpreted as a suppression of instanton effects by the light quark modes ||. Of the three 
quantities discussed in this section, x shows by far the most significant unquenching features. 
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Figure 2: Topological susceptibility, x, versus (effective) quark mass. All quantities have been scaled 
by the distance scale r^. 
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4 Spectral Functions 



The 2-point function, G^i) can De spectrally decomposed 

POO 

G 2 (t) = / K(t,uj)p(uj)duj, 
Jo 

where K is the Kernel (it is e~ ujt in our case - see Eq(|l])) and p(u>) is the Spectral Function. p(co) 
contains information about the ground and excited states and the continuum for the particular 
channel of interest. 

In H the "Maximum Entropy Method" which is based on Bayesian statistics was employed to 
extract p(co) directly from the lattice G2(t) data. This is an exciting development which we plan 
to implement since a huge wealth of information becomes available once p{oj) is calculated. This is 
especially true in QCD in the high temperature phase where many channels no longer have bound 
states, and the traditional approach (detailed in section 2) is simply not applicable. 



5 Summary 

The UKQCD collaboration have made important progress in lattice QCD by implementing a 
"matched" trajectory approach where the lattice spacing is kept fixed, and thus the unquench- 
ing, finite-size and discretisation effects are dis-entangled. A summary is given here of the results 
for the static quark potential, hadronic spectrum and topological susceptibility - see [||, ||] for full 
details. A new and very promising method which enables the direct extraction of the spectral 
function from lattice data is discussed. 
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